arXiv:1504.06788v2 [gr-qc] 14 Aug 2015 


X-ray spectropolarimetric measurements of the Kerr metric 


Dan Liu, Zilong Li, Yifan Cheng, and Cosimo Bambj^] 

Center for Field Theory and Particle Physics and Department of Physics, Fudan University, 200433 Shanghai, China 

(Dated: August 17, 2015) 

It is thought that the spacetime geometry around black hole candidates is described by the Kerr 
solution, but an observational confirmation is still missing. Today, the continuum-fitting method 
and the analysis of the iron K a line cannot unambiguously test the Kerr paradigm because of the 
degeneracy among the parameters of the system, in the sense that it is impossible with current 
X-ray data to distinguish a Kerr black hole from a non-Kerr object with different values of the 
model parameters. In this paper, we study the possibility of testing the Kerr nature of black hole 
candidates with X-ray spectropolarimetric measurements. Within our simplified model that does 
not include the effect of returning radiation, we find that it is impossible to test the Kerr metric and 
the problem is still the strong correlation between the spin and possible deviations from the Kerr 
geometry. Moreover, the correlation is very similar to that of other two techniques, which makes the 
combination of different measurements not very helpful. Nevertheless, our results cannot be taken 
as conclusive and, in order to arrive at a final answer, the effect of returning radiation should be 
properly taken into account. 
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I. INTRODUCTION 

General relativity is today the best framework for the 
description of the gravitational force and the geometri¬ 
cal structure of the spacetime around massive bodies. In 
the past 60 years, the theory has passed a number of ex¬ 
perimental tests and its predictions have been verified in 
the weak gravitational fields, mainly with precise exper¬ 
iments in the Solar System and accurate radio observa¬ 
tions of binary pulsars [T. Now the interest is to check 
the validity of the theory in more extreme environments. 
One of the most fascinating predictions of general relativ¬ 
ity is the existence of black holes (BHs). In 4-dimensional 
general relativity, an uncharged BH is described by the 
Kerr solution and it is completely characterized by only 
two parameters, namely its mass M and its spin param¬ 
eter a* = a/M = J/M 2 , where J is the BH spin angular 
momentum. Kerr BHs are thought to be the final stage 
of any heavy star after it exhausts all its nuclear fuel [2]. 
Initial deviations from the Kerr solution can indeed be 
quickly radiated away through the emission of gravita¬ 
tional waves [3| • Any non-vanishing electric charge would 
be soon neutralized because of the highly ionized host en¬ 
vironment of these objects [4]. Deviations from the Kerr 
metric produced by the presence of accretion disks are 
normally completely negligible, as the disk mass is typi¬ 
cally many orders of magnitude smaller than the mass of 
the central body [5]. 

Astrophysical BH candidates are dark and compact 
objects that can only be interpreted as Kerr BHs in the 
framework of conventional physics and they can be some¬ 
thing else only in the presence of new physics. For in¬ 
stance, a compact object in an X-ray binary is classified 
as a BH candidate if its mass exceeds 3 M©, because the 
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latter is the maximum mass for a neutron star for any rea¬ 
sonable matter equation of state |6]. Present observations 
cannot confirm that the spacetime geometry around BH 
candidates is really described by the Kerr metric [7 . For 
the time being, there are only two relatively robust tech¬ 
niques to study the nature of BH candidates, namely the 
analysis of the thermal spectrum of thin accretion disks 
(continuum-fitting method) [8j and the iron Ka line [9]. 
Both techniques have been developed to estimate the spin 
parameter under the assumption of the Kerr background 
and, more recently, they have been extended to test the 
nature of BH candidates [T014T2] . 

The typical problem to verify the Kerr BH hypothesis 
is the degeneracy among the parameters of these systems, 
and in particular among the estimate of the spin, possi¬ 
ble deviations from the Kerr geometry, and the inclina¬ 
tion angle of the disk. The disk’s thermal spectrum has a 
simple shape and therefore it is fundamentally impossible 
to distinguish the effect of the spin from non-Kerr met¬ 
ric elements m- The iron line has a more complicated 
structure and, in the presence of the correct astrophysical 
model, it would be possible to distinguish Kerr and non- 
Kerr BHs with high quality data. However, this is impos¬ 
sible with current X-ray facilities, even in the case of very 
good observations [T4] . Moreover, some kinds of defor¬ 
mations are definitively more difficult to constrain than 
others PI- With current data of the iron line, we can 
only rule out some BH alternatives without horizon, like 
some boson stars and some traversable wormholes [US] , 
because their iron line profile would have qualitatively 
different features. We note that the combination of the 
measurement of the continuum and of the iron line of the 
same object can unlikely break the degeneracy, because 
both techniques are mainly sensitive to the position of the 
inner edge of the disk m, which is normally supposed 
to be at the radius of the innermost stable circular orbit 
(ISCO) and it is thus only determined by the background 
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metricQ The studies of quasi-periodic oscillations and of 
jet power are potentially other techniques to probe the 
metric around BH candidates |18], but the exact mech¬ 
anism responsible for these phenomena is not yet well 
understood and therefore they cannot yet be used to test 
fundamental physics. 

The study of the polarization of the thermal radia¬ 
tion of thin accretion disks may be a new technique to 
test the Kerr metric of stellar-mass BH candidates in the 
near futur^J Such a radiation is initially unpolarized, 
but it gets polarized at the level of a few percent due to 
Thomson scattering of X-ray radiation off free electrons 
in the disk’s atmosphere. In the Kerr metric, the degree 
and the angle of polarization depend on the BH spin and 
the inclination angle of the disk with respect to the line 
of sight of the observer m- Assuming the Kerr back¬ 
ground, spectropolarimetric observations could provide 
and estimate of these two parameters [20] 21 • While 
some polarimetric missions have been cancelled, similar 
measurements will be hopefully possible in the near fu¬ 
ture with the Chinese X-ray Timing Polarimetric (XTP) 
satellite [2]2], which may be launched in 2020, and the 
European X-ray Imaging Polarimetry Explorer (XIPE) 
mission [23], which was recently approved for initial de¬ 
sign work. 

X-ray spectropolarimetric measurements can be poten¬ 
tially used to test the Kerr metric. Such a possibility has 
been already explored in Ref. [24], where is was found 
that these measurements are mainly sensitive to the po¬ 
sition of the ISCO. In the present paper, we want to per¬ 
form a more detailed analysis on the correlation among 
the parameters of the system and on future detection ca¬ 
pabilities. As an explorative study in this direction, we 
make some simplifications, and in particular we neglect 
the returning radiation, namely the radiation emitted by 
the disk that comes back to the disk because of the strong 
light bending in the vicinity of the compact object. First, 
we confirm the strong correlation between the spin and 
possible deviations from the Kerr geometry. Second, the 
correlation is very similar to that found in the estimates 
with the continuum-fitting method and the iron line anal¬ 
ysis, which suggests that the possible combination of the 
three measurements is not promising to break the degen¬ 
eracy. Third, considering some plausible future measure¬ 
ments of the polarization of the continuum, we find that 
constraints with this method are not better than those 
from the thermal spectrum of thin disks and surely worse 
than those from the iron line. Even in the case of objects 


1 We note that it is often assumed that the emission of the plunging 
gas inside the ISCO can be completely neglected. 

2 Here we are interested in stellar-mass BH candidates because the 
disk temperature scale as M -1 / 1 2 * 4 and in the case of a compact 
object of 10 Mq the peak of the spectrum is around 1 keV. For 
supermassive BH candidates of millions or billions Solar masses, 

the thermal spectrum of a thin disk falls in the UV/optical band, 

where dust absorption makes an accurate measurement impossi¬ 
ble. 


that look like fast-rotating Kerr BH observed from quite 
a large inclination angle, it is impossible to exclude sig¬ 
nificant deviations from the Kerr metric. Fourth, if we 
do not assume a priori the Kerr background, we lose also 
the capability of obtaining stringent constraints on the 
inclination angle of the disk. 

We note that our results are not conclusive, because we 
have neglected the effect of returning radiation. The lat¬ 
ter depends on the light bending in the vicinity of the BH 
and it is possible that its signature improves the measure¬ 
ment of the Kerr metric, even because it is not strictly 
related to the position of the ISCO radius. As clearly 
shown in Ref. [21], the signature imprinted by the return¬ 
ing radiation is weak below 10 keV for a Schwarzschild 
BH (roughly speaking, when the ISCO radius is not too 
close to the compact object). However, it becomes promi¬ 
nent for a fast-rotating BH from energies above a few 
keV, which is an energy range covered by the proposed 
polarimetric missions that should do measurements at 
^1-10 keV. The effect will be included in a future work, 
since it cannot be treated as in the Kerr background and 
the calculations become very time consuming. 

The content of this paper is as follows. In the next 
section, we briefly review our theoretical model for the 
description of the accretion disk and the spacetime met¬ 
ric of a BH candidate. In Section m we describe our 
calculations of the spectrum of the polarization degree 
and the polarization angle. In Section [IV] we study the 
degeneracy among the parameters of our model, namely 
the spin parameter, the deformation parameter, and the 
inclination angle of the disk with respect to the line of 
sight of the distant observer. Summary and conclusions 
are reported in Section 0 In the appendix, we provide 
some details about the calculations of the polarization of 
the radiation from the disk, since we cannot use the usual 
approach adopted in the Kerr metric exploiting the nice 
properties of the Kerr geometry. Our approach is more 
general, in the sense it can be applied to any stationary, 
axisymmetric, and asymptotically flat spacetime, but it 
is inevitably more time consuming. Throughout the pa¬ 
per, we use units in which Gn = c = 1 . 


II. THEORETICAL MODEL 
A. Accretion disk 

We describe the accretion disk with the Novikov- 
Thorne model [25] . which is the standard set-up for ge¬ 
ometrically thin and optically thick accretion disks in 
stationary and axisymmetric spacetime [8 . The disk is 
assumed to be on the equatorial plane, the particles of the 
gas follow nearly geodesic circular orbits, and the inner 
edge of the disk is at the ISCO radius. See e.g. the last 
paper in |8] for the validity of these assumptions. From 
the conservations laws for rest-mass, energy, and angular 
momentum, we can derive the time-averaged structure 
of the disk. The disk is in thermal equilibrium and the 
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emission at any radius is like that of a blackbody, so that 
we can define an effective temperature as a function of 
the radial coordinate, T e ff(r). Since the temperature in 
the inner part of the accretion disk around a stellar-mass 
BH candidate can be very high, up to ~ 10 7 K, electron 
scattering in the atmosphere makes the spectrum devi¬ 
ate from the original Novikov-Thorne prediction. This 
can be taken into account by introducing the color factor 
(or hardening factor) f c [26]. Since we are interested in 
stellar-mass BH candidates with an accretion luminosity 
of about 10% the Eddington luminosity of the object, 
throughout the paper we use f c = 1.6. The color tem¬ 
perature is defined as T c = f c T e $ and the local specific 
intensity of the radiation emitted by the disk is 

%) = 1 b(t c )t, (i) 

J c 

where v e is the photon frequency in the rest frame of the 
gas, B(T C ) is the blackbody function for the temperature 
T = T c , and T is a function of the angle between the 
direction of the propagation of the photon and the normal 
to the disk surface. 

Thermal radiation is initially unpolarized. However, 
because of the Thomson scattering of photons off free 
electrons in the dense atmosphere of the disk, the radia¬ 
tion becomes partially polarized. With reference to the 
rest frame of the gas, the degree of polarization depends 
on the angle d between the normal to the disk surface and 
the direction of propagation of the X-ray photon, ranging 
from 0 {d = 0 °, photon direction parallel to the normal 
to the disk) to about 12% (d = 90°, photon direction 
orthogonal to the normal to the disk) [27]. The polariza¬ 
tion degree as a function of d (which is measured in the 
rest frame of the gas) is shown in Fig. [l] The orientation 
of the polarization vector is instead parallel to the disk 
plane and orthogonal to the direction of propagation of 
the photon. The same scattering in the disk atmosphere 
causes a limb-darkened emission and the correct value for 
T can be found in the table in m- More details can be 
found in MEO. 

In the present paper, we adopt two main simplifica¬ 
tions in the calculation of the polarization. We do not 
take into account the effect of the returning radiation, 
namely the effect of the photons that are emitted by the 
disk and, because of the strong gravitational field near 
the compact object, return to the disk, increasing its tem¬ 
perature and local specific intensity [ 2 T] . Such an effect 
is clearly more important in the high energy part of the 
spectrum, because high energy photons are produced at 
small radii, and for fast-rotating BHs, because the inner 
edge of the disk is closer to the compact object. In the 
latter case, the effect can indeed be important just above 
a few keV (for more details, see Ref. ( 21 |) and therefore 
our results cannot be taken as conclusive. 

The second ingredient that is neglected in our calcu¬ 
lation is the photon absorption, which tends to destroy 
the polarization and it is more important for low energy 
photons, so at larger radii ED]. This effect is expected 



FIG. 1. Polarization degree as a function of the angle d be¬ 
tween the normal to the disk surface and the direction of 
propagation of the photon. 


to be important below ~ 0.5 keV and therefore should 
not be relevant for the first generation of polarimetric 
detectors. 


B. Background metric 


In order to test the Kerr metric around BH candidates, 
it is necessary to constrain possible deviations from the 
Kerr geometry. Indeed, it is not enough that observa¬ 
tional data nicely fit a Kerr model, because a non-Kerr 
object may look like a Kerr BH with different values of 
the model parameters. As example, we can consider the 
approach used in Solar System experiments. In this case, 
we want to check the Schwarzschild solution in the weak 
field limit. In the Parametrized Post-Newtonian (PPN) 
formalism |1 , we write the most general static and spher¬ 
ically symmetric metric as an expansion in M/r 


, 2 , 2 M 2 M 2 

ds 2 = - 1 -+ 


dt 2 


+ ^1 + ( dx 2 + dy 2 + dz 2 ) , ( 2 ) 


where /? and 7 are two coefficients in the expansion and 
they parametrize our ignorance. In general relativity, 
the only spherically symmetric vacuum solution is the 
Schwarzschild metric and, when cast in the form above, 
we have (3 = 7 = 1. To test the Schwarzschild metric, 
we employ the line element in Eq. ©> and we determine 
/3 and 7 from observations. Current data require that (3 
and 7 are 1 at the level of 10 “ 5 — 10 -4 and this confirms 
the Schwarzschild solution with this precision [l]. 

The same strategy can be employed to test the Kerr 
metric. At present, there is not a satisfactory approach 
like the PPN formalism: since we want to probe the 
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spacetime close to the compact object, we cannot use 
an expansion in M/r and it is thus difficult to take into 
account any kind of deviations from the Kerr solution. In 
this paper, we adopt the Johansenn-Psaltis metric [28] , 


which is a quite popular metric to test the Kerr space- 
time. In Boyer-Lindquist coordinates, the line element 
reads [28] 


ds 2 = - 1 - 


2 Mr 


, _. _ 9 4 aMr sin 2 6 , 

(1 + h)dt 2 --- (1 


h ) dtd<p ■ 


E(l + fo) 


-dr 2 


+Xd<9 2 + 


r 2 + a 2 + 


2a 2 Mr sin 2 0\ . 2 a 2 (X + 2Mr) sin 4 0 
-—- sm 2 0 -\ ---——- h 


A + a 2 hsin 2 0 
4^ 

# 2 ■ 


(3) 


where X = r 2 + a 2 cos 2 6 and A = r 2 — 2 Mr + a 2 . In its 
simplest version, h is given by 


e 3 M 3 r 

'*—25- 


(4) 


63 is the “deformation parameter” and it is used to quan¬ 
tify possible deviations from the Kerr geometry. The 
compact object is more prolate (oblate) than a Kerr BH 
for 63 > 0 (63 < 0 ); when 63 = 0 , we exactly recover the 
Kerr solution. 


III. NUMERICAL SIMULATIONS 

We use the code described in Ref. m, which has been 
extended to compute the degree and the angle of po¬ 
larization of the thermal spectrum of a thin disk. We 
remind the reader that we are performing these calcula¬ 
tions in a non-Kerr background, and therefore we cannot 
adopt the usual approach used in a Kerr code, in which 
one exploits the fact that the Kerr solution is a Petrov 
type D spacetime. As discussed in m , the calculations 
of the photon propagation backward in time from the ob¬ 
server’s plane to the point of the photon emission in the 
accretion disk is done by solving the geodesic equations. 
Now we also need to compute the polarization degree S 
and the polarization angle ^ for any photon on the ob¬ 
server’s plane. The polarization degree is a scalar and it 
only depends on the angle $ between the normal to the 
disk surface and the direction of propagation of the X-ray 
photon and it does not require special prescription. For 
the polarization angle, we need to parallel transport the 
polarization vector along the photon geodesic; we can¬ 
not exploit the Walker-Penrose theorem valid for Petrov 
type D spacetimes, and therefore we need to solve the ba¬ 
sic equations for parallel transport. At any point of the 
grid of the observer’s plane, we simultaneously compute 
the photon trajectory (backward in time) and the prop¬ 
agation of an auxiliary vector. When the photon reaches 
the disk, we evaluate the polarization degree and the an¬ 
gle difference between the propagated auxiliary vector 
and the polarization vector at the emission point. Since 
the angle between the two vectors is conserved along the 
geodesic, we can immediately determine the angle of the 


polarization vector on the plane of the distant observer. 
The details are given in the appendix at the end of this 
paper. 

Once we know the polarization degree and angle at 
each point in the grid, we need to integrate over the ob¬ 
server’s plane to get the spectrum of S and ^ (h * 1 this 
paper we use the notation of Ref. [20]). In terms of the 
Stokes parameters /, Q, U, and V [27] . for each point on 
the image we have 


Q + iU = 6Ie 2i * , (5) 


where V = 0 because the radiation is linearly polar¬ 
ized. The radiation field is decomposed into a completely 
polarized component I p = SI and an unpolarized one 

I u = (1 — 5)1. At the point of the detection, we have 

(Qobs) 4" i{U Q bs) — A 0 f (Qobs ffi iU Q bs) d£l obs 
^^obs J 

= a?S— / 9 3 S e I e e 2i ^dn obs , (6) 
A^obs J 

where (•) indicates the average over the image, the 
subindices “obs” and “e” refer, respectively, to quanti¬ 
ties measured in the rest-frame of the observer and of the 
emitter, AU 0 b s is the total solid angle subtended by the 
disk on the sky, and the redshift factor g = E 0 \> s /E e en¬ 
ters from the conservation of the quantity I/E 3 along the 
photon path, namely J 0 bs/^obs = ^e/^e- We note that, 
in general, the component of the radiation that is ini¬ 
tially completely polarized is detected on the observer’s 
plane as partially polarized, because different points of 
the image have photons with different V’obs- The total 
intensity at the detection point is 

(A)bs) = TV [9 7edU 0 bs 

^^obs J 

= <4“bs) + <4 P bs>- (7) 

The observed averaged polarization degree is m 


{Sobs) 


V^Qobs ) 2 + (t/obs ) 2 

(r> bs) 


(8) 


and the observed averaged polarization angle is deter- 
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FIG. 2. Polarization degree (left panel) and polarization angle (right panel) as a function of the photon energy for a 
Schwarzschild BH (solid lines) and a Kerr BH with spin parameter a/M — 0.9 (dashed lines) and a viewing angle i = 45°, 60°, 
and 75°. See the text for more details. 
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FIG. 3. Apparent image of the accretion disk around a Kerr 
BH with spin parameter a/M — 0.99 and observed from a 
viewing angle i = 75°. The contour map shows the total 
intensity of the radiation (logarithmic scale). The segments 
report the properties of the polarization of the radiation: the 
length of every segment is proportional to the polarization 
degree, while the orientation of every segment corresponds to 
the orientation of the polarization vector. See the text for 
more details. 


mined from the following two relations [2CI 


sin (2 (V’obs)) = 


(C obs ) 

y(Qobs) 2 + (Cobs) 2 ’ 


(Cobs) 

yiCobs ) 2 + ' 


( 9 ) 


With the above machinery, we can compute the spec¬ 
trum of the polarization degree and of the polarization 
angle for a specific model. In our non-Kerr model, there 
are six basic parameters, namely the mass M, the dis¬ 
tance d, the inclination angle of the disk with respect to 
the line of sight of the observer i, the mass accretion rate 
M, the spin parameter a*, the deformation parameter 
63. In what follows, we always assume M = 10 M 0 and 
an accretion luminosity to Eddington luminosity ratio of 
0.1. The latter sets the mass accretion rate via L = r]M, 
where 77 is the radiative efficiency in the Novikov-Thorne 
model (see the discussion in Section 3 in Ref. [13]). We 
have checked the results of our code in the Kerr met¬ 
ric with those reported in Ref. [21] . Fig. [2] shows the 
spectrum of the polarization degree (left panel) and the 
spectrum of the polarization angle (right panel) for a 
Schwarzschild BH (solid lines) and a Kerr BH with spin 
parameter a* = 0.9 (dashed lines) when the viewing an¬ 
gle is i = 45°, 60°, and 75° (from bottom to top in the 
left panel, from top to bottom in the right panel). These 
plots can be compared with those in Fig. 2 in Ref. [2T] 
and the agreement is good. The apparent image of the 
accretion disk around a Kerr BH with spin parameter 
a* = 0.99 is shown in Fig. [3J The contour map shows 
the relative intensity of the total (namely polarized and 
unpolarized) radiation, I Q b s /^obs,max (logarithmic scale). 
The black segments show the polarization of the radia¬ 
tion: the length of the segment is proportional to the 
polarization degree, while its orientation corresponds to 
that of the polarization vector. Our Fig. [ 3 ] can be com¬ 
pared with Fig. 1 in Ref. [21] and the result is definitively 
similar. 


COS (2('0obs)) — 


( 10 ) 
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FIG. 4. Ay 2 contour lines on the plane spin parameter a* and deformation parameter 63 from X-ray spectropolarimetric 
measurements. The reference model is a Kerr BH with spin parameter a* = 0.9 and observed from a viewing angle i = 70°. In 
the left panel, all models use an inclination angle of 70°. In the right panel, the inclination angle is free in the fit. See the text 
for more details. 




FIG. 5. As in Fig. [4] on the plane inclination angle i and deformation parameter £ 3 . In the left panel, all models use a spin 
parameter a* = 0.9. In the right panel, the spin is free in the fit. See the text for more details. 


IV. TESTS OF THE KERR METRIC 

In this section, we want to figure out how X-ray spec- 
tropolarimetric measurements of the thermal spectrum 
of a thin disk (the so-called continuum component) can 


constrain the Kerr metric around a BH candidate, and, in 
particular, the parameter degeneracy. For a systematic 
study, we start considering a model with specific values 
of the spin parameter, deformation parameter, and incli¬ 
nation angle (reference model). This is compared with a 
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FIG. 6 . Ay 2 contour lines on the plane spin parameter a* and deformation parameter 63 from X-ray spectropolarimetric 
measurements. The reference model is a non-Kerr BH with spin parameter a* = 0.65, deformation parameter e 3 = 4, and 
observed from a viewing angle i — 70°. In the left panel, all models use an inclination angle of 70°. In the right panel, the 
inclination angle is free in the fit. See the text for more details. 



FIG. 7. As in Fig. [ 6 ] on the plane inclination angle i and deformation parameter 63 . In the left panel, all models use a spin 
parameter a* = 0.65. In the right panel, the spin is free in the fit. See the text for more details. 


set of Kerr and non-Kerr models with different values of 
a*, €3, and i (template models) with the use of x 2 - As 


done in Ref. [21 , we define x 2 as 


X 2 (a*, 63 , i) = X 


k =1 L 


(Qk-Qf ) 2 + (u k -ufy 


A Q‘l 


At4 2 


,(n) 



























where the summation is performed over n sampling en¬ 
ergies Ek, Qk and Uk are the Stokes parameters of the 
template spectrum with parameters a*, 63, and i in the 
energy bin [Ek , Ek + A E ], while Q )^ ef and U £ ef are the 
Stokes parameters of the reference spectrum in the en¬ 
ergy bin [Ek , E k + A E]. The errors A Qk and A Uk are 
assumed to be 

IJ ref 

A Q k = A U k = S min lfJ , (12) 

where S m [ n is the minimum polarization sensitivity at the 
peak of the spectrum, J£ ef is the total (polarized and un¬ 
polarized) intensity of the reference spectrum in the en¬ 
ergy bin and k is the total intensity of the reference 
spectrum at the peak. In this way we see how an obser¬ 
vation can constrain the parameters of a putative source 
with the same values of the spin, the deformation, and 
the inclination angle as the reference model. 

In this work, we assume that the energy range of the 
detector is 0.5-10 keV, the energy resolution is A E = 
0.1 keV, and that S m i n = 0.003. The results of our sim¬ 
ulations are reported in Figs. [4][7] These plots show the 
contour levels Ay 2 = y 2 — Xmin = 1, 4, 9, and 16, where 
here Ay 2 = y 2 because Xmin = 0? since the minimum 
corresponds to the reference model without noise (the in¬ 
troduction of the noise would not change our results, of 
course). In the case of one degree of freedom, the 68.3% 
C.L. designed as 1-standard deviation limit corresponds 
to Ay 2 = 1 [29] • In the case of three degrees of freedom, 
we have Ay 2 = 3.53, 8.03, and 14.16, respectively for 
68.3%, 95.4%, and 99.7% (1-, 2-, 3-standard deviation 
limit). 

As the first case, as reference model we consider a Kerr 
BH with spin parameter a* = 0.9 and observed from an 
inclination angle i = 70°. We note that a similar ob¬ 
ject would be a good source for our purpose, because 
the value of both the spin and the inclination angle are 
quite high and this should maximize relativistic effects, 
making the difference in the spacetime geometry more 
evident. Fig. [4] and Fig. [5] show, respectively, the con¬ 
straints on the (a*, €3) and (i, 63) planes. In the left pan¬ 
els, the third parameter is fixed, while in the right panels 
it is free in the fit and we minimize y 2 with respect to 
it. From Fig. [4| we see that there is a strong correla¬ 
tion between the estimate of a* and 63 and that it is 
very similar to that found in the case of the continuum¬ 
fitting method and the iron line analysis (see in particular 
Refs. [13, IT?]). The possibility of an independent accu¬ 
rate measurement of i is not very helpful to break the 
degeneracy between a* and 63, and indeed the difference 
between the left and right plots is small. This conclusion 
is confirmed by Fig. [5] An independent good measure¬ 
ment of a* would be enough to constrain €3 (left panel), 
but without it large deviations from the Kerr geometry 
cannot be excluded (right panel). 

Fig.© and Fig. [7] show the contour levels for a hypo¬ 
thetical non-Kerr BH with a * = 0.65 and 63 = 4. The 
inclination angle of the reference model is still i = 70°. 


These plots confirm the conclusion reported in the pre¬ 
vious paragraph. We should note that in Fig. [6] the al¬ 
lowed region is actually much longer than that appar¬ 
ently shown and that it extends to the Kerr models with 
63 = 0. The allowed region is very narrow for a fixed 
63, and therefore the numerical simulations miss it, but 
this can be understood by noting that both in the left 
and right panels there are some allowed “islands” near 
a* = 1 and 63 = 0. Once again, the Kerr metric could 
be tested with an independent estimate of a* (Fig. [7| left 
panel), while it is not possible to do it if a* is free in the 
fit (Fig.0 right panel). 


V. SUMMARY AND CONCLUSIONS 

Current observations cannot unambiguously confirm 
whether astrophysical BH candidates are the Kerr BHs 
or general relativity because of a degeneracy among the 
parameters of the system: a non-Kerr compact object 
may mimic a Kerr BH with a different spin parameter 
and observed from a slightly different viewing angle. In 
the case of the shape of the thermal spectrum of thin 
disks, such a degeneracy is at the fundamental level and 
more accurate measurements cannot break it [TO] . Con¬ 
cerning the iron line profile, this technique is potentially 
more powerful: with the available X-ray data, the intrin¬ 
sic Poisson noise of the source does not allow to break 
the degeneracy between the spin and possible deviations 
from the Kerr geometry, but future high-quality data can 
do it m, even if not for any kind of deformations m- 

The radiation of the thermal spectrum of thin accre¬ 
tion disks around BH candidates is inevitably polarized 
due to Thomson scattering of X-ray photons off free elec¬ 
trons in the dense atmosphere above the disk. The de¬ 
gree of polarization depends on the photon direction with 
respect to the normal of the disk surface, while the po¬ 
larization vector is initially parallel to the plane of the 
disk. However, the strong gravitational field in the vicin¬ 
ity of BH candidates strongly affects the degree and the 
angle of polarization detected by an observer far from 
the object and polarization measurements can provide 
information about the metric. In this paper, we have 
studied the possibilities offered by X-ray spectropolari- 
metric measurements to test the Kerr geometry around 
BH candidates. Our targets are stellar-mass BH candi¬ 
dates in the high/soft state. This kind of measurements 
may be possible in the near future with the XTP [22] and 
the XIPE [23] missions. 

As the first step in our program, we have considered 
a simple theoretical model, in which we neglect the ef¬ 
fects of returning radiation and photon absorption by the 
disk’s atmosphere, two ingredients that will be included 
in future developments. We have considered an hypo¬ 
thetical detector with an energy resolution of 0.1 keV in 
the 0.5-10 keV band and a detection sensitivity of the 
polarization degree at the peak of the spectrum of 0.003. 
As reference models, we have studied a Kerr BH with 
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spin a* = 0.9 and a non-Kerr BH with spin a* = 0.65 
and deformation parameter €3 = 4. Both the reference 
models have an inclination angle i = 70°. Our results 
can be summarized as follows: 

1. We confirm the conclusions of Ref. |24], namely 
that X-ray spectropolarimetric measurements can¬ 
not unambiguously test the Kerr metric because 
there is a strong correlation between the estimate 
of the spin and possible deviations from the Kerr 
solution. 

2. The correlation between these two parameters and 
the estimate of the inclination angle of the disk is 
instead modest, with the result that an indepen¬ 
dent accurate measurement of i does not help to 
break the degeneracy between the spin and the de¬ 
formation parameter. 

3. The correlation between the estimates of the spin 
and the deformation parameter in X-ray spectropo¬ 
larimetric measurements is very similar to those al¬ 
ready found in the case of the continuum-fitting 
method and of the analysis of the iron Ka line. 
This suggests that the possible combinations of the 
three measurements is not promising to break the 
degeneracy. 

4. In the case of a plausible future polarization mea¬ 
surement, it seems that the capability to test the 
Kerr metric with this technique is roughly com¬ 
parable to that offered by the continuum-fitting 


method and definitively worse than high quality 
data of the iron line. 

5. If we relax the Kerr BH hypothesis, constraints on 
the disk’s inclination angle with X-ray spectropo¬ 
larimetric measurements become weaker. 

The characteristic of our hypothetical detector are prob¬ 
ably a little bit too optimistic for a first-generation of 
X-ray spectropolarimetric detectors, and therefore such 
a technique may not compete with the available ones in 
the near future. However, we must note that our results 
cannot be taken as conclusive. As shown in Ref. EQ, 
the scattering of returning radiation can imprint its sig¬ 
nature just above a few keV in the case of fast-rotating 
Kerr BHs. Since the effect is caused by the light bend¬ 
ing in the vicinity of the object, it may provide further 
constraints and/or an independent information with re¬ 
spect to the thermal spectrum or the iron line. A final 
answer on the possibility of testing the Kerr metric with 
X-ray polarimetric data thus requires to properly include 
the effect of returning radiation, and this is left to future 
work. 
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Appendix A: Polarization calculations 


We consider an observer at a distance d from the BH and with an inclination angle i. His/her image plane is provided 
with a system of Cartesian coordinates (A, Y, Z\ as shown in Fig. [8] Another system of Cartesian coordinates (x, y , z) 
is centered at the BH with the orientation reported in Fig. [8j The two Cartesian coordinates are related by 

x = dsini — Y cosi + Zsini , y = X, z = D cosi + Ysinz + Zcosi . (Al) 

Far from the BH, the Boyer-Lindquist spatial coordinates reduce to the usual spherical coordinates and they are 
related to (x, y , z) by 


r = \/ x 2 + y 2 + z 2 , 


6 = arccos 



(j) = arctan 



(A2) 


1. Initial conditions 

The image of the accretion disk is provided by the photons hitting the image plane of the distant observer with 
3-momentum perpendicular to the plane. In our numerical calculations, we start from photons on the image plane 
and we follow their trajectory and the evolution of their polarization vector till they hit the accretion disk. 
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FIG. 8. Set-up of the system. The Cartesian coordinates (x,y,z) are centered at the BH. The image plane of the distant 
observer is located at the distant d from the BH, with an inclination angle i, and it is provided with a system of Cartesian 
coordinates (X, Y, Z). In the calculations, a photon starts from the image plane with initial position (Xo,Yo,0), initial 3- 
momentum ko = —koZ perpendicular to the image plane, and auxiliary 3-vector fo = X perpendicular to its 3-momentum and 
parallel to the X-axis. See the text for more details. 


Let us consider a photon at the position (X 0 , Yq, 0), with 3-momentum k 0 = — k$Z perpendicular to the image 
plane, and an auxiliary 3-vector fo = X. The initial conditions for the photon position are 


*0 = 0, r 0 = JX 2 + Y 2 + d 2 , 9 0 = 


Yq sin i + d cos i 

-F- 


, (j) o = arc tan 


X 0 


d sin i — Yq cos i ’ 


(A3) 


and for the photon 4-momentum are 


d 

nr 


k o = -|ko|, Ko 


fcn = 


cos i — (lo sin i + d cos i) -4 
— ^g- l 

\/Xq + (d sini — Fq cosi) 2 


mj| • H ) 


W = 


Xq sin i 


Xq + (d sin i — Yq cos i) 2 


(A4) 


and fcg = \J{Kf + Tq (^) 2 + Tq sin 2 #o(&o) 2 f°U° ws from the condition g^k^k 1 ' = 0 with the metric tensor of a flat 
space-time. The initial conditions for the auxiliary vector are 


f 0 = 0, f5 = —, /? = 

Tq 


d cos i + Yq sin i 


Xq 


yjX q + (d sin i — Yq cos i) 2 r ° 


2 ’ 


f<p _ 

Jo — 


d sin i — Yq cos i 


Xq + {d sin i — Yq cos i) 


■ \ 2 ' 


(A5) 


2. Propagation in non-Kerr spacetimes 


The photon trajectory is found by solving the geodesic equations 


d 2 x* dx"dx° 

d\ 2 ua dX dX 


(A6) 


and the auxiliary vector is parallel transported along the photon geodesic 


d f* fV dx^_ 

dX vaI dX 


(A7) 


In order to have numerical errors under control, at any step we check that the following relations are satisfied: 

/4/4 = i, Fk^ = o. 
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3. Coordinate transformations at the emission point 

When the photon hits the disk, we have to determine: i) the angle id between the normal to the disk and the 
direction of propagation of the photon as measured in the rest frame of the accreting gas, ii) the orientation of the 
auxiliary vector as measured in the rest frame of the accreting gas. The angle id is necessary to obtain the polarization 
degree. The orientation of the auxiliary vector is compared to that of the polarization vector expected from the 
emerging radiation. Since the angle between the two vectors is conserved along the photon path, it is the same at 
the emission and at the detection points, and therefore we can immediately determine the angle of the polarization 
vector measured far from the BH. 

Our spacetimes are stationary, axisymmetric, and asymptotically flat. We write the line element in the following 
form 

ds 2 = — e 2u dt 2 + e 2fJj dr 2 + e 2X d0 2 + e 2cr (dcf) — ujdt) 2 . (A8) 

The covariant basis vectors in the locally non-rotating frame is 

E^= (e",0,0,0), £«= (0,e^,0,0), E™ = (o, 0, e A , o) , E™ = ( - , 0, 0, ) . (A9) 

The gas 4-velocity in the locally non-rotating frame is = E^u 1 ' and therefore the velocity of the gas with respect 
to the locally non-rotating frame is 

U W 

v=— v = (n-c j) e a ~ u , (A10) 

uW 

where Q = /u 1 is the angular velocity of the gas as measured by the distant observer. The covariant basis vectors in 

the rest-frame of the gas are related to those in the locally non-rotating frame by the following Lorentz transformation 

E® = 7 {e® - vE<E) , E<p = £« , E& = E™ , = 7 (_ vE W + _ ( A11) 

where 7 = (l — v 2 ) is the Lorentz factor. Eventually, we can write the photon 4-momentum and the polarization 
4-vector in the rest-frame of the gas 

jfe( fl ) “ , / (o) = Eft F *. (A12) 


4. Calculation of d and ^ 


The polarization degree S depends on the angle id between the normal to the disk surface and the direction of 
propagation of the photon. Once we know $, we can find the corresponding S from the table in Ref. m- Since the 
emission is on the equatorial plane, we have 


cos'd = 






y(fcw) 2 + (fcw) 2 + 




id = arccos 


m 


(A13) 


For the calculation of the polarization angle, we proceed in the following way. The polarization vector of the 
radiation emerging from the disk, say E a \ is oriented perpendicular to the propagation direction of the photons and 
parallel to the disk (vanishing d-component). We choose the gauge in which the ^-component of the polarization 
vector vanishes and, from the normalization condition = 1 and the orthogonality with the photon momentum 

^(a) = 0, we find 

^ / 1 \ k 

(A14) 


h^ r) = 


= - 


( U r + [k^y 


(k( r ^ + (kwy 


Before comparing with , we need to use the same gauge for the two vectors and we thus remove the ^-component 
in f( a \ namely 

/ ,( °) = / (a) - ak < a ) , 


/ (a) 


(A15) 


r 


where a = /k^\ We define the angle £ as 


and the polarization angle detected by the distant ob- 


£ = 


arccos ( f r h r + 


if f > 0 


arccos (— f r h r — if 7 < 0 


(16) 
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server is 




{I 


i-t 


if^<§ 

if e > f 


(17) 
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